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We present a neutron diffraction study of the magnetic structure of single crystal PrC>2 under ap- 
plied fields of 0-6 T. As the field is increased, changes are observed in the magnetic Bragg intensities. 
These changes are found to be irreversible when the field is reduced, but the original intensities can 
be recovered by heating to T > 1 22 K, then re-cooling in zero field. The antiferromagnetic ordering 
temperature Tn = 13.5 K and the magnetic periodicity are unaffected by the applied field. We also 
report measurements of the magnetic susceptibility of single crystal PrC>2 under applied fields of 0- 
7T. These show strong anisotropy, as well as an anomaly at T = 122±2K which coincides with the 
temperature Td = 120 ±2 K at which a structural distortion occurs. For fields applied along the [100] 
direction the susceptibility increases irreversibly with field in the temperature range Tn < T < Td. 
However, for fields along [110] the susceptibility is independent of field in this range. We propose 
structural domain alignment, which strongly influences the formation of magnetic domains below 
Tn, as the mechanism behind these changes. 

PACS numbers: 61.12.Ld, 75.25.+Z, 75.30.Cr, 75.30.Kz 



I. INTRODUCTION 

In recent years there has been strong interest in Jahn- 
Teller and orbital phenomena in compounds containing 
localized 4/ and 5/ electrons. Among the simplest of 
these are the fluorite-structure actinide dioxides UO2 and 
Np02 , which display complex ordered phases at low tem- 
peratures involving coupled electric and magnetic multi- 
poles as well as (in the case of UO2) a lattice distortionJa 2 " 

Unusual magnetic effects have also been observed in 
the lanthanide dioxide Pr02, which is isostructural with 
UO2 and Np02 at room temperature and exhibits anti- 
ferromagnetic ordering below Tn = 13.5 K. Some years 
ago, Pr02 was found to have an anomalously small or- 
dered moment in the antiferromagnetic phased More re- 
cently we discovered a broad continuum in the magnetic 
excitation spectrum probed by neutron inelastic scatter- 
ing, which we ascribed to Jahn- Teller fluctuations in- 
volving the orbitally degenerate 4/ ground state and 
dynamic distortions of the latticed In a separate neu- 
tron diffraction experiment!* we found evidence that the 
antiferromagnetic structure contains a component with 
twice the periodicity of the accepted magnetic structure, 
and we have recently reported further neutron diffraction 
studies* which reveal an internal distortion of the fluorite 
structure below Td = 120 ± 2K and a related distortion 
of the antiferromagnetic structure below Tn- These dis- 
tortions result in a doubling of both the crystallographic 
and magnetic unit cells along one crystal axis. The mag- 
netic structure is found to consist of two components: a 
primary component whose unit cell is the same as that 
of the undistorted crystal structure (referred to hereafter 
as the "type-I component" ) , and a secondary component 
with a smaller ordered moment, whose unit cell is the 



same as that of the distorted structure (referred to here- 
after as the "doubled component"). 

In this paper we present neutron diffraction studies 
of the magnetic structure and measurements of the mag- 
netic susceptibility of Pr02 under applied magnetic fields 
of 0-7 T. The neutron diffraction experiments were de- 
signed to influence the alignment of symmetry-equivalent 
magnetic domains, in order to determine whether the 
type-I component of the magnetic structure was of the 
multi-q type. These experiments produced unexpected 
results which suggest that the alignment of the magnetic 
domains is strongly influenced by a field-induced align- 
ment of structural domains in the distorted phase below 



II. MAGNETIC STRUCTURE UNDER 
APPLIED FIELD 

As mentioned in the introduction, the magnetic struc- 
ture of Pr02 contains two components, the dominant 
one being the type-I antiferromagnetic structure which 
has a magnetic unit cell of the same dimensions as the 
unit cell of the undistorted fluorite crystallographic struc- 
ture. The secondary (doubled) component of the mag- 
netic structure is described in detail elsewhere»S There 
are three possible type-I spin arrangements consistent 
with existing neutron diffraction data. These are the 
transverse multi-q structures shown in Fig. ^ Longi- 
tudinal structures are ruled out by comparison of their 
magnetic structure factors with the measured Bragg in- 
tensities. Neutron diffraction measurements in zero field 
cannot distinguish between these three structures be- 
cause they give rise to identical magnetic structure fac- 
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tors when averaged over symmetry-equivalent magnetic 
domains. One of the aims of the experiments described 
here was to influence the populations of the domains by 
cooling through Tn in a magnetic field, and hence to dis- 
tinguish between the three structures. 
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FIG. 1: Transverse multi-q magnetic structures for Pr02- 
The spheres are the Pr ions (the O ions are not shown in this 
diagram) . 
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FIG. 2: Field dependence of the (011) magnetic Bragg reflec- 
tion before (open circles) and after (closed circles) the irre- 
versible increase in zero-field intensity. The field was applied 
parallel to [Oil] and the temperature was held constant at 
T = 1.55 K during both field sweeps. 



A. Experimental Details 

Two neutron diffraction experiments were performed 
using the same crystal in different orientations with re- 
spect to the magnetic field. The first experiment was 
performed on the E4 double-axis single crystal diffrac- 
tometer at the Berlin Neutron Scattering Centre at the 
Hahn-Meitner Institute. A flat pyrolytic graphite (002) 
monochromator was used in combination with a graphite 
filter. 40' collimators were placed before and after the 
monochromator, but there was no collimation between 
the sample and the detector. A square aperture of 10 mm 
x 10 mm was placed before the detector. The incident 
neutron wavelength was 2.44 A. 

The second experiment was performed on the 
D10 single-crystal diffractometer at the Institut Laue- 
Langevin. This was operated in double-axis mode with a 
position-sensitive detector. A vertically curved pyrolytic 
graphite (002) monochromator was used in combination 
with a graphite filter. No collimators were used, but a 
circular aperture of diameter 6 mm was placed in the in- 
cident beam before the sample, and a square aperture 
of 12 mm x 12 mm was placed before the detector. The 
incident neutron wavelength was 2.356(4)1. 

For both experiments we used a single crystal sample of 
Pr02 which was prepared some time ago by a hydrother- 
mal procedure^ and had a mass of ~ 1 mg. For the first 
experiment the crystal was mounted inside a vertical-field 
superconducting cryomagnet with a temperature range of 
2-300 K and a magnetic field range of 0-5 T. It was glued 
onto a thin aluminium pin such that the [011] direction 
was vertical and hence parallel to the applied field. For 
the second experiment the crystal was mounted inside a 
similar cryomagnet with a slightly larger field range of 
0-6 T. It remained mounted on the aluminium pin, but 
was re-aligned using an attachment that held it at 45° so 
that the [001] direction was vertical. 



B. Measurements with H || [011] 

We first describe measurements performed on the E4 
diffractometer. Five magnetic reflections of the type-I 
component of the magnetic structure were accessible in 
the horizontal scattering plane of the crystal ([011] verti- 
cal). These were (100), (011), (211), (122) and (300). At 
the start of the experiment the crystal was cooled to a 
temperature T = 1.55 K with zero applied field, and these 
five magnetic reflections plus several structural Bragg re- 
flections were measured by crystal rotation (w-scan). 

After this we slowly increased the field from H = 0T 
to H = 5 T in 0.5 T steps (while remaining at T = 1.55 K) 
and monitored the count at peak centre of the (011) re- 
flection. The intensity decreased continuously with field, 
reducing almost to zero at H = 5 T. As we did not warm 
above Tn during the field sweep the disappearance of the 
(011) reflection suggested a spin reorientation transition 
rather than a change in the populations of the magnetic 
domains. To check that we could regain the original in- 
tensity, we immediately warmed the sample to T = 20 K 
under H = 5T, removed the field, cooled the sample 
back to T = 1.55 K in zero field, and re-measured the 
five type-I magnetic Bragg peaks. Surprisingly, the orig- 
inal intensities were not recovered, w-scans revealed that 
the intensities of all five peaks had increased by ~ 100% 
relative to the initial zero- field measurement. The inten- 
sities of the (200) and (022) structural Bragg peaks were 
checked and found to be the same as they had been before 
application of the field. We then repeated the field scan 
from H = 0Tto J ff = 5TatT = 1.55 K while counting 
at the centre of the (011) peak. The intensity was again 
found to reduce to zero (see Fig.|2). By monitoring the 
zero-field intensity of the (100) peak as a function of tem- 
perature, we determined that the Neel temperature had 
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not changed, remaining at T N = 13.5 K. 

To illustrate the increase in intensities of the five type-I 
magnetic Bragg peaks we compare in Fig.|3]the (100) and 
(Oil) peaks, measured in approximately zero field, before 
and after application and removal of the 5 T field asso- 
ciated with the above temperature cycling. We should 
mention that when the peaks were remeasured after the 
irreversible increase in intensities the magnetic field was 
actually set to H = 0.5 T, rather than H = 0T, but sub- 
sequent field scans showed almost no difference between 
0T and 0.5 T (see Fig.©. 



„ 1000 

"o 900 

1 800 
■ 700 

| 600 

| 500 

2 400 
I 300 

200 



Before * 


(100) - 


(H = T) 


• After 


H//[0T1] 


- (H-0.5T) % 

4 4 

,•«.«•»♦ 


' T = 1.55K - 



-78.0 -77.5 -77.0 -76.5 -76.0 
i» (degrees) 



550 

500 
» 450 
ii 400 

| 350 
| 300 
«T 250 

1 200 



( 



t 



(011) 



Before 
(H = 0T) 

After H // [01 1 ] 

(H = 0.5T) j T=1.55K 



! ( 



1 * * i 



62.5 -162.0 -161.5 -161.0 -160.5 
w (degrees) 



FIG. 3: w-scans of the (100) and (011) magnetic Bragg re- 
flections at T — 1.55 K, before (open circles) and after (closed 
circles) application and removal, above Tn, of a 5T magnetic 
field parallel to the [Oil] direction. Detector counts are nor- 
malised to a fixed incident beam monitor count of 4 x 10 J , 
corresponding to a counting time of ~ 4 minutes. 



In an attempt to recover the original intensities of the 
five peaks the sample was heated to T — 35 K and cooled 
back down to 1.55 K in zero field. However, the origi- 
nal intensities were not recovered. The sample was then 
heated to T = 123 K and cooled back down in zero field. 
This time the original intensities were recovered. 

As mentioned in the introduction, PrC>2 undergoes a 
structural distortion at T D = 120 ± 2 K, which is accom- 
panied by an anomaly in the magnetic susceptibility. 6 
The fact that we were able to recover the original in- 
tensities of the magnetic Bragg reflections by heating 
to T > T D suggests that the PrC>2 ground state is 
metastable in the distorted phase. From now on we will 
refer to the Bragg intensities obtained by cooling through 
T D in zero field as the "original intensities" . 

The intensities of three of the reflections, (100), (011) 
and (211), were also measured by w-scan at H = 5T 
(after cooling through T N in the 5 T field). The (100) and 
(211) reflections were found to have increased by ~ 100% 
from their original intensities at H = 0T, whereas the 
(011) reflection was found to have decreased almost to 
the level of the background as reported above (see Fig. 
I2J). The (100) and (011) peaks are shown in Fig. H for 
comparison with those shown in Fig. [21 

Since the (100) peak did not disappear at H = 5T, 
we decided to check its field dependence. Starting from 
the original intensity, we found that when the field was 
applied at T = 1.55 K there was little change in intensity 
up to H — 5T. However, when the crystal was warmed 
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FIG. 4: oj-scans of the (100) and (011) magnetic Bragg reflec- 
tions (normalised to a fixed incident beam monitor count), 
measured at H — 5 T after cooling through Tn in the field 
(H || [Oil]). 



above Tn and cooled to 1.55 K each time the field was 
incremented we observed a smooth increase in intensity 
with field. Furthermore, the increase did not become 
irreversible until the applied field was larger than H = 
3.5 T. Once the 100% increase had been achieved, the 
intensity of the (100) reflection was unaffected by the 
applied field, whether the crystal was cooled through Tn 
in the field or not. The same was found to be true of the 
(211) and (300) reflections. However, the (122) reflection 
behaved more like the (011) reflection, decreasing to zero 
intensity at H = 5 T when cooled through Tn in the 
applied field. 

The main findings from the measurements described 
so far with H || [011] are (i) the disappearance of certain 
magnetic Bragg reflections on application of a field H ~ 
5T, (ii) an irreversible 100% increase in intensity of all 
five peaks when the field is removed above T N and the 
sample re-cooled to 1.55 K, (hi) recovery of the original 
intensities after heating above the structural distortion 
temperature T D = 120 K and re-cooling in zero field. We 
also note that the irreversible increase in intensities can 
be generated by applying a field H ~ 5T at 20 K, then 
cooling through T N either before or after removing the 
field. The application of fields up to 5 T at T = 1.55 K 
does not cause the zero-field intensities to increase. 



C. Measurements with H || [001] 

The experiment performed on the D10 diffractometer 
was a repeat of the study performed on E4, but this 
time with the field applied parallel to the [001] direction. 
Since the field was constrained to be vertical in both ex- 
periments the crystal orientation had to be changed for 
the D10 experiment. This altered the scattering plane 
and hence the range of accessible magnetic reflections. 
The type-I reflections accessible with our chosen neutron 
wavelength were (010), (110), (120), (030), (130) and 
(230). A number of half- integer reflections corresponding 
to the doubled component of the magnetic structure were 
also accessible, but only the (^10) and (|10) reflections 
had sufficient intensity to be measurable. 
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We measured these eight peaks by w-scan at T = 2.4 K 
in zero applied field. We then warmed to 18 K, applied 
a field of 5.7 T, and re-cooled through Tn in the ap- 
plied field. Before re-measuring the peaks we warmed 
to 18 K and removed the field before cooling back to T = 
2.4 K (this second warming/cooling cycle was performed 
to reproduce the conditions under which the irreversible 
increase in intensities was first observed in the E4 ex- 
periment, i.e. field removed above Tjj). The repeat w- 
scans revealed that the (110) and (130) peaks reduced by 
~ 80%, while the other type-I peaks increased by ~ 40%. 
The (A 10) and (§10) peaks increased by - 20% and 30% 
respectively. Figure shows the intensities of the (010) 
and (110) peaks before and after application and removal 
of the 5.7 T field associated with the above temperature 
cycling. 
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FIG. 6: Field dependence of the (010) and (110) magnetic 
Bragg reflections at peak centre before (open symbols) and 
after (closed symbols) the irreversible change in peak inten- 
sities. The magnetic field was applied parallel to the [001] 
direction, and the crystal was warmed above Tn and cooled 
back to T — 2.4 K each time the field was changed. Note that 
the monitor count for this data is four times that of the data 
displayed in Figure |S] 
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FIG. 5: The intensities of the (010) and (110) magnetic Bragg 
reflections at T = 2.4 K, measured at H — 0T before (open 
circles) and after (closed circles) application and removal, 
above Tn, of a 5.7 T magnetic field parallel to the [001] di- 
rection. A monitor count of 1 x 10 6 corresponds to a counting 
time of ~ 2 minutes. 

We then measured the (010), (110) and (230) peaks 
at T = 2.4 K, H > 5T after cooling through T N in the 
field. Their intensities remained approximately the same 
as those observed at zero field after application and re- 
moval of the 5.7 T field. 

Figure [fj] shows a plot of the count at peak centre for 
the (010) and (110) reflections as a function of applied 
field before and after the change in intensities. Each time 
the field was changed the crystal was warmed above Tn 
and cooled back to T = 2.4 K. Once the maximum field 
of 5.7 T had been applied, further changes in the applied 
field had little effect on the peak intensities. 

In agreement with expectations following the E4 exper- 
iment, we found that application and removal of a large 
magnetic field (5.7 T) above Tn, followed by zero-field 
cooling caused an irreversible change in the intensities 
of the magnetic Bragg peaks. We also found that the 
original intensities could only be regained by heating the 
crystal to T > Td and re-cooling through Tb in zero field. 
The intensities of the (020) and (220) structural Bragg 
reflections were checked at T = 2.4 K, H = T before and 
after application of the 5.7 T field. Again, these did not 
change. There was no change in the Neel temperature 
either, this remaining at Tn = 13.4 K for both integer 
and half-integer magnetic reflections. 



The most striking result from the measurements de- 
scribed above is that the changes in the peak intensi- 
ties are different from those observed during the E4 ex- 
periment. In the E4 experiment (when H was applied 
along the [011] direction) all the peak intensities rose 
by ~ 100% after application and removal of a 5 T field, 
whereas in the D10 experiment (when H was applied 
along the [001] direction) the relative peak intensities 
changed: some peaks increased by ~ 40% while others 
decreased by ~ 80%. 

D. Data analysis 

We obtained the integrated intensity of each magnetic 
Bragg reflection measured in the experiments described 
above by fitting a Gaussian profile to each peak, calcu- 
lating its area and correcting for the geometrical Lorentz 
factor L = sin 2$, where 9 is the Bragg angle. Tables H] 
and UTI list the integrated, corrected intensities and com- 
pare them with the square of the expected magnetic 
structure factor (|Tm(Q)| 2 ) for a type-I antiferromag- 
netic structure under ambient conditions, i.e. where all 
the magnetic domains are equally populated. The mag- 
netic structure factors were calculated using the formula 

<|Tm(Q)| 2 ) = £((<W - Q a Qp) FS(Q)<(Q)) , 

a/3 

where the summation indices a and /3 run over the carte- 
sian co-ordinates x, y and z, 8 a p is the Kronecker delta, 
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Q a is the a-component of the unit scattering vector and 
() denotes an average over all symmetry-equivalent mag- 
netic domains. F^(Q) is given by 



Fu(Q) 



/(Q)Vi?e iQ ^ 



3 



-W 3 (Q,T) 



(2) 



where the summation index j runs over all the magnetic 
atoms in the magnetic unit cell, /t" is the a-component 
of a unit vector in the direction of the magnetic moment 
of the jth magnetic atom, Yj is the position of the jth 
magnetic atom within the magnetic unit cell, /(Q) is the 
magnetic form factor of the Pr 4+ ion and e~ w ^Q-^ is 
the Debye- Waller factor (we set this equal to 1, since the 
measurements were made at low temperatures). 

The integrated intensities of the magnetic Bragg re- 
flections are proportional to (|-Fm(Q)| 2 ) and the volume 
of the sample. For ease of comparison the intensities of 
the observed magnetic reflections listed in tables [I] and 
UTI have been normalised to the original intensity of the 
(100) reflection. Similarly, the (|_Fm(Q)| 2 ) have been nor- 
malised to the (|F M (Q)| 2 ) of the (100) reflection. 



Reflection 


H || [011] 


Domain- 
averaged 

(I^m(Q)I 2 ) 


Before 


After 


H = 


H = 0.5 T 


H = 5T 


(100) 


1.00 ± 0.05 


1.87 ± 0.05 


1.93 ± 0.06 


1.00 


(011) 


0.56 ± 0.04 


1.00 ± 0.06 


0.09 ± 0.05 


0.48 


(211) 


0.63 ± 0.07 


1.34 ± 0.09 


1.18 ± 0.09 


0.65 


(122) 


0.47 ± 0.07 


0.91 ± 0.08 


*0.51 ± 0.26 


0.38 


(300) 


0.56 ± 0.07 


1.17 ± 0.09 


*1.67 ± 0.31 


0.69 



TABLE I: Comparison between integrated intensities of mag- 
netic peaks (corrected for the Lorentz factor) before and after 
application of a 5T magnetic field along [Oil]. The inten- 
sities have been normalised to the original intensity of the 
(100) reflection. The (|F M (Q)| 2 ) have been normalised to the 
(|Jm(Q)| 2 ) of the (100) reflection. The values marked with 
a * have been calculated from measurements of the count at 
peak centre, since no cj-scans were made at these positions. 

For H || [011] the relative intensities of the mag- 
netic Bragg reflections agree with the domain-averaged 
(|-Fm(Q)| 2 ) both before and after the irreversible increase 
in intensities. However, for H || [001] the relative peak in- 
tensities only agree with the domain-averaged (|-Fm(Q)| 2 ) 
before the irreversible change. 



E. Summary 

To summarise the results of the neutron diffraction 
studies on single crystal PrC>2, we have observed an irre- 



Reflection 


H || [001] 


Domain- 
averaged 

(\Fm{Q)\ 2 ) 


Before 


After 


H = 


H = 


H w 5T 


(010) 


1.00 ± 0.09 


1.43 ± 0.10 


1.38 ± 0.10 


1.00 


(110) 


0.50 ± 0.06 


0.10 ± 0.03 


0.09 ± 0.05 


0.48 


(120) 


0.53 ± 0.11 


0.76 ± 0.10 




0.49 


(030) 


0.57 ± 0.08 


0.79 ± 0.08 




0.69 


(130) 


0.26 ± 0.06 


0.07 ± 0.02 




0.33 


(230) 


0.47 ± 0.07 


0.65 ± 0.08 


0.76 ± 0.16 


0.49 


(|io) 


0.08 ± 0.02 


0.10 ± 0.02 






(fio) 


0.10 ± 0.02 


0.13 ± 0.03 







TABLE II: Comparison between integrated intensities of mag- 
netic peaks (corrected for the Lorentz factor) before and after 
application of a 5.7 T magnetic field along [001]. The inten- 
sities have been normalised to the original intensity of the 
(010) reflection. The (|-Fm(Q)| 2 ) have been normalised to the 
(|F M (Q)| 2 ) of the (010) reflection. 



versible change in the magnetic Bragg intensities follow- 
ing the application and removal of a large magnetic field. 
The field required to induce this change is between 3.5 T 
and 5 T when H || [011] and < 5.7 T when H || [001]. To 
recover the original intensities the crystal must be heated 
to T > Tb and re-cooled through Tb in zero field. The ex- 
perimental results suggest that the field must be applied 
above Xn to induce the irreversible change in intensities. 
When the field is applied parallel to the [011] direction, 
certain peaks increase by ~ 100% while others disappear. 
After removal of the field, all peaks are found to have in- 
creased by ~ 100% relative to their original intensities. 
When the field is applied parallel to [001] some of the 
type-I peaks increase by ~ 40% while others decrease 
by ~ 80%. The half-integer peaks increase by 20-30%. 
These changes are preserved when the field is removed. 
The intensities of the structural Bragg reflections are un- 
affected by fields applied along either the [011] direction 
or the [001] direction, and Tn remains unchanged. 



III. INFLUENCE OF MAGNETIC FIELD ON 
DOMAIN POPULATIONS 

We now attempt to interpret the observations sum- 
marised above by considering the possible influence of an 
applied magnetic field on the populations of structural 
and magnetic domains. 

One of the original aims of the neutron diffraction ex- 
periments described above was to influence the popula- 
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tions of magnetic domains by cooling through Tn in a 
magnetic field. By applying the field along certain sym- 
metry directions we hoped to break the symmetry of the 
crystal in such a way that certain domains would be en- 
ergetically favoured. In the single-q and doublc-q mag- 
netic structures some domains have zero magnetic struc- 
ture factors. Therefore, altering the populations of the 
domains can alter the magnetic Bragg intensities. How- 
ever, in the triple-q structure all domains have identi- 
cal magnetic structure factors, so in this case the Bragg 
intensities will be unaltered by changes in the domain 
populations. 

The fact that we have observed changes in the inten- 
sities of the magnetic Bragg peaks after application of 
a magnetic field along two different symmetry directions 
rules out the triple-q structure. We can also rule out the 
double-q structure by calculating the single-domain mag- 
netic structure factors for the observed reflections. We 
find no single domain whose squared magnetic structure 
factor is more than 50% larger than its value when av- 
eraged over all symmetry-equivalent domains. Thus the 
double-q structure cannot account for the 100% increase 
in intensities observed when the field is applied along 
[Oil]. However, to prove that the structure is single-q 
and that the increase in intensities observed is due to 
preferential population of certain domains we must de- 
termine whether the increased intensities are consistent 
with the magnetic structure factors of the domains we 
expect to be favoured. 

The single-q transverse^ structure has six symmetry- 
equivalent domains, as the ordering wavevector can be 
along [100], [010] or [001], and the spins can point along 
either of the two directions mutually perpendicular to 
the ordering vector. Tables ITTT1 and llVl show the single- 
domain |i 7 lvi(Q)| 2 calculated in these six domains for the 
magnetic Bragg peaks observed in the E4 and D10 exper- 
iments respectively. We use to denote a domain with 
ordering vector along the a-direction and spins along the 
/^-direction. 

If we assume that the exchange interactions between 
the individual atomic spins are stronger than their inter- 
actions with the applied field, then for an antiferromag- 
netic structure we expect the domains with the lowest 
magnetic energy to be those whose spins lie perpendicu- 
lar to the applied field. We therefore expect such domains 
to be favoured over those whose spins have components 
parallel and antiparallcl to the field. If we use this as- 
sumption to predict which magnetic domains should be 
favoured in PrC>2, we expect domains and D° to be 
favoured equally for H || [011]. However, the magnetic 
structure factors in these domains are zero for all the 
peaks measured in the E4 experiment, in contradiction 
to the observed irreversible increase in intensities. Sim- 
ilarly, for H || [001] we would expect domains D b a , D%, 
D° and D h c to be favoured equally. However, if we aver- 
age the magnetic structure factors for these domains we 
find that they do not agree with the observed intensity 
changes. 



Reflection 


Single-domain li^Q)! 2 


Domain-averaged 


D b a 


Dt 


Dt 


D% 




D b c 


(100) 


3.00 


3.00 


0.00 


0.00 


0.00 


0.00 


1.00 


(011) 


1.42 


1.42 


0.00 


0.00 


0.00 


0.00 


0.48 


(211) 


1.96 


1.96 


0.00 


0.00 


0.00 


0.00 


0.65 


(122) 


1.14 


1.14 


0.00 


0.00 


0.00 


0.00 


0.38 


(300) 


2.06 


2.06 


0.00 


0.00 


0.00 


0.00 


0.69 



TABLE III: Comparison between single-domain and domain- 
averaged |.Fm(Q)| 2 for magnetic Bragg reflections observed 
in the E4 experiment (where the magnetic field was applied 
along [011]). The single-domain |-Fm(Q)| 2 have been nor- 
malised to the domain-averaged (|-Fm(Q)| 2 ) of the (100) re- 
flection. Z)f denotes a domain with ordering vector along the 
a-direction and spins along the /3-direction. 



Reflection 


Single-domain -F M (Q) 2 


Domain-averaged 

(I^m(Q)| 2 ) 


D b 




D a b 


Dt 


£>c 


D\ 


(010) 


0.00 


0.00 


3.00 


3.00 


0.00 


0.00 


1.00 


(110) 


0.00 


0.00 


0.00 


0.00 


1.42 


1.42 


0.48 


(120) 


0.49 


2.46 


0.00 


0.00 


0.00 


0.00 


0.49 


(030) 


0.00 


0.00 


2.06 


2.06 


0.00 


0.00 


0.69 


(130) 


0.00 


0.00 


0.00 


0.00 


1.77 


0.20 


0.33 


(230) 


0.00 


0.00 


1.20 


1.74 


0.00 


0.00 


0.49 



TABLE IV: Comparison between single-domain and domain- 
averaged |_Fm(Q)| 2 for magnetic Bragg reflections observed 
in the D10 experiment (where the magnetic field was ap- 
plied along [001]). The single-domain |i<M(Q)| 2 have been 
normalised to the domain-averag cd (|F M (Q)| 2 > of the (010) 
reflection. denotes a domain with ordering vector along 
the a-direction and spins along the /3-direction. 



The above observations suggest that our assumptions 
concerning the dominant spin interactions are incorrect. 
By inspection of table UTTl we find better agreement be- 
tween the structure factors and increased intensities if do- 
mains D b a and D c a are favoured for H || [011] . Rather than 
having spins perpendicular to the applied field, these do- 
mains have ordering vectors perpendicular to the applied 
field. If these domains are favoured, the intensities of all 
the observed peaks should increase by a factor of 3. The 
observed increase is a factor of ~ 2, but it is possible 
that this is due to incomplete depopulation of the dis- 
favoured domains. Similarly, by inspection of table IIVI 
we find better agreement between the structure factors 
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and observed intensities if domains D° a , and 
are favoured for H || [001]. Again, these domains have 
ordering vectors (rather than spins) perpendicular to the 
applied field. If we average their squared magnetic struc- 
ture factors we find that the (110) and (130) peaks are 
expected to disappear, while the rest are expected to in- 
crease by 50%. The observed decrease of the (110) and 
(130) intensities by ~ 80% and the increase of the other 
peaks by ~ 40% are consistent with a small portion of 
the crystal remaining populated with the disfavoured do- 
mains. 

In the D10 experiment we also observed the (^10) and 
(flOj reflections to increase by 20-30% after applica- 
tion and removal of a large field parallel to [001]. These 
reflections arise from the doubled component of the mag- 
netic structured If we calculate the single-domain mag- 
netic structure factors of the (|l0) and (§10) reflections, 
we find that they are non-zero only in domains where 
the unit cell is doubled along the a-direction. Since 
the field was applied along [001], i.e. the c-direction in 
the D10 experiment, we suggest that the applied field 
favours domains of the doubled component of the mag- 
netic structure whose unit cells are doubled along the 
a- and 6-directions, i.e. directions perpendicular to H. 
The squared magnetic structure factors of the QlO) and 
(| 10) reflections averaged over these domains are 50% 
larger than the structure factors averaged over all do- 
mains, in rough agreement with the observed increase in 
intensities of the peaks. 

We have shown that magnetic domain alignment can 
account for the observed changes in magnetic Bragg in- 
tensities following application and removal of a large 
magnetic field. However, we have not yet discussed why 
the changes are irreversible below Tr> We propose that 
this is because certain domains of the distorted crystal 
structure are energetically preferred in an applied field. 
In other words, the application of a magnetic field at con- 
stant temperature below Tb causes a change in the struc- 
tural domain population. Intuitively, it feels somewhat 
implausible that a 5 T magnetic field could cause a signif- 
icant repopulation of structural domains at low temper- 
atures. However, we must remember that the structural 
distortion is purely internal, so does not produce any 
macroscopic strain. All that is required is that the oxy- 
gen atoms move the small distance between equivalent 
shifted sites within the unit cell. Since it is likely that 
the cooperative Jahn- Teller distortion at Td is accom- 
panied by ordering of the Pr orbitals, we propose that 
the mechanism by which the field influences the struc- 
tural domains is a coupling between the applied field and 
the Pr orbital magnetic moments. Our analysis indicates 
that the favoured structural domains are those whose 
unit cells are doubled along directions perpendicular to 
the applied field. The existence of a doubled component 
of the magnetic structure indicates a strong coupling be- 
tween the lattice and the magnetic ordering, possibly via 
the Pr spin-orbit interaction. We therefore propose that 
the formation of magnetic domains on cooling through 



Tn is influenced predominantly by the underlying align- 
ment of the structural domains. 

From our analysis of the favoured type-I magnetic do- 
mains it appears that coupling between the type-I compo- 
nent and the doubled component of the magnetic struc- 
ture causes the ordering vector of the type-I component 
to lie parallel to the direction along which the unit cell 
is doubled. The strong influence of the structural do- 
main alignment on the magnetic ordering of both the 
type-I component and the doubled component would ac- 
count for the irreversibility of the changes in the magnetic 
Bragg intensities below Td- 

Finally we recall that in the E4 experiment the (011) 
and (122) reflections disappeared under the influence of 
a 5T magnetic field along [011]. This occurred both 
when the field was applied at constant T = 1.55 K and 
when the crystal was cooled through Tn in the field. The 
reflections reappeared when the field was removed. To 
explain these observations we propose that a reversible 
spin-reorientation transition occurs in the type-I compo- 
nent between H = T and H = 5 T for H || [Oil]. This 
causes the spins to rotate parallel and antiparallel to the 
field as shown in Fig. \7\ We arrived at this spin con- 
figuration by trying various possibilities until we found 
one whose squared single-domain magnetic structure fac- 
tors were in agreement with the relative intensities of 
the magnetic peaks observed at H = 5 T (see Table [Vj . 
The configuration appears counterintuitive because an 
antiferromagnetic structure would normally minimise its 
energy by rotating its spins perpendicular to the field. 
However, our experimental results indicate a strong in- 
teraction between the lattice and the magnetism which 
could alter this simple picture. 




A 

H II [Oil] 



FIG. 7: Reorientated spin configuration of the type-I compo- 
nent of the magnetic structure at H = 5 T. 



IV. SUSCEPTIBILITY MEASUREMENTS 

We now present measurements of the magnetic suscep- 
tibility of single crystal PrC>2. These measurements were 
made using a commercial SQUID magnetometer with a 
vertical field range of 0-7 T and a temperature range of 
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Reflection 


Intensity at H — 5 T 


Smgle-domam _F M (Q) 


(100) 


1.93 ± 0.06 


3.00 


(Oil) 


0.09 ± 0.05 


0.00 


(211) 


1.18 ± 0.09 


1.57 


(122) 


*0.51 ± 0.26 


0.23 


(300) 


*1.67 ± 0.31 


2.05 



TABLE V: Comparison between integrated intensities of mag- 
netic peaks (corrected for the Lorentz factor) measured un- 
der a field of 5T parallel to [Oil] and single-domain |_Fm(Q)| 2 
for the structure shown in Fig. [7| The intensities have been 
normalised to the original intensity of the (100) reflection. 
The |Jm(Q)| 2 have been normalised to the domain-averaged 
(|-Pm(Q)| 2 ) of the (100) reflection for the type-I magnetic 
structure. The values marked with a * have been calculated 
from measurements of the count at peak centre, since no uo- 
scans were made at these positions. 



2-300 K. The susceptibility was measured as a function 
of temperature and field, with the field applied along two 
different symmetry directions. The aim was to probe the 
bulk magnetisation of the crystal as a function of mag- 
netic field to observe the effects of domain alignment and 
spin reorientation. 



A. Experimental Details 



Figures Efa) and |8£b) show the temperature depen- 
dence of the magnetic susceptibility in the range T = 
2-300 K when a field of H = 1 T is applied parallel to 
the [110] direction and when a field of H = 0.3 T is ap- 
plied parallel to the [001] direction. Figures[Sfc) and[Sf d) 
show the same data as in Figs. Ufa) and|SJb) respectively, 
but depict the inverse susceptibility as a function of tem- 
perature, which makes the anomaly at Td = 120 K more 
visible. The overall shape of the susceptibility trace is 
very similar for H || [110] and H || [001]. However, the 
anomaly at Td is more pronounced for H || [110]. This 
can be seen in the plots of both M/H vs T and (M/7T) -1 
vs T. 
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We used two single crystals of PrC>2 for the suscep- 
tibility measurements. Both were taken from the same 
batch as the sample used in the neutron diffraction ex- 
periments, and each had a mass of < 1 mg. The crystals 
were mounted in plastic sample holders with the [110] 
and [001] directions vertical (hence parallel to the ap- 
plied field). Their orientations were estimated to be ac- 
curate to within 5%. All measurements were made using 
the reciprocating sample option (RSO) which causes the 
sample holder to oscillate in the vertical direction. 



B. Measurements and results 

Here we present measurements of the magnetic sus- 
ceptibility performed with the applied field parallel to 
the [110] and [001] directions. Results obtained for the 
two different field directions are compared and discussed. 

It should be noted that the susceptibility measure- 
ments presented in this section are not normalized to the 
mass of the crystals because the crystals were too small 
to be weighed accurately. Hence, the susceptibility data 
taken in different field directions should not be compared 
quantitatively. 



FIG. 8: Temperature dependence of the magnetic suscepti- 
blity of PrC>2 when the magnetic field is applied parallel to 
the [110] and [001] directions. Plots (c) and (d) show the 
inverse susceptibilities derived from the data in (a) and (b) 
respectively. 



To investigate the field-dependence of the susceptibil- 
ity, we made a series of measurements between 2 K and 
20 K at different applied fields. The field was increased 
and decreased in steps from ff = 0Ttoif = 7T. For each 
measurement we applied the field above Tn (at 20 K), 
then measured the susceptibility as a function of temper- 
ature while cooling through Tn. The data are shown in 
Fig. El with field increasing in[3fa) and^Jc) and decreas- 
ing in|5Jb) and^d). For both field directions the suscep- 
tibility below Tn increases smoothly with field through 
the pale grey region (see^a) andHJc)). Plots |!Jc) and 
EJd) are identical, showing that, for H || [110] the increase 
below Tn is completely reversible and the susceptibility 
above Tn is independent of field. For H || [001], a change 
occurs at H = IT. Below this field the susceptibility 
above Tn is independent of field, but above H = 1 T 
the whole trace increases irreversibly through the dark 
grey region shown in[!|{a). When the field is removed, 
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the whole trace decreases through the smaller dark grey 
region shown in[§fb). As the field is decreased below 
H = 1 T the susceptibility above Tn becomes indepen- 
dent of field and the trace below Tn decreases through 
the pale grey region. 
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FIG. 10: Change in the field-dependence of the susceptibility 
with temperature. These plots are taken from the same data 
set as those in Fig.|^] Each time the applied field was changed, 
the crystal was warmed to T = 20 K, then cooled through Tn 
in the new applied field, (a) and (b) show M/H for H || [001] 
below and above Tn respectively, (c) and (d) show M/H for 
H || [110] below and above T N . 



Temperature (K) 



Temperature (K) 



FIG. 9: Field-induced changes in the temperature- 
dependence of the susceptibility. The field is increased and 
decreased in steps, with each new field applied at T — 20 K. 
The susceptibility is measured as a function of T as the crys- 
tal is cooled through Tn. (a) and (b) show M/H for H || 
[001]. The field increases in (a) and decreases in (b). (c) and 
(d) show M/H for H || [110]. The field increases in (c) and 
decreases in (d). 



It should be noted that repetitions of the susceptibil- 
ity measurements at low applied fields revealed small, 
random shifts whenever the applied field was changed 
between measurements. This was attributed to fluctu- 
ations in the magnetisation of the sample surroundings 
and trapped flux in the superconducting magnet, lead- 
ing to fluctuations in the applied field of 0.001-0.02 T. 
However, at higher fields the fluctuations became negli- 
gible compared to the total field, so the susceptibility was 
unaffected. Unfortunately it was impossible to measure 
the magnetisation of the sample surroundings, and it was 
therefore impossible to measure the exact applied field. 
This meant that the magnitude of the susceptibility at 
applied fields of H < 0.2 T had an uncertainty of at least 
± 10%. 

The plots in Fig. ^| are taken from the same data 
set as those in Fig. [5J but this time the susceptibility is 
plotted as a function of applied field instead of tempera- 
ture. Figures and llOf c) show the field dependence 
of the susceptibility below Tn, whereas Figs. HuT b') and 
UUd) show the field dependence above T N . For H || [001] 
the susceptibility below Tn increases rapidly from H = 

0- 1 T, then increases slowly and irreversibly from H = 

1- 7T (Fig. IToT all. The susceptibility above Tn changes 



little between H = T and H = 1 T, but increases irre- 
versibly between H = 3 T and H = 7T fFig-HOTb^. For 
H || [110] the susceptibility below Tn increases rapidly 
and reversibly between H = 0T and H = 7T, but is 
almost independent of field above Tn- 

In the measurements described above, the field was 
always applied above Tsr- However, we also wished to 
investigate the field dependence of the susceptibility be- 
low Tn without cooling through Tn for each change in 
field. For this measurement the magnetic field was swept 
while keeping the temperature constant at 2K. Figure 
II II shows the results obtained. For H || [110] the trace is 
very similar to that obtained when the crystal was cooled 
through Tn in the field. However, for H || [001] there are 
two differences. Firstly, the irreversible increase occurs 
at a higher field. Secondly, on removal of the field at T 
= 2K, the susceptibility does not return to its original 
zero field value. This contrasts with the data shown in 
Fig. llOf a - ). where removal of the field above Tn reduced 
the susceptibility below its original zero field value. 

To see if an irreversible increase in the susceptibility 
occurred at higher temperatures we measured the field 
dependence at T = 60 K, with H || [001]. An irreversible 
increase did occur (see Fig. I12f a)). but of smaller mag- 
nitude than that observed at T = 20 K. To find the field 
required to make the increase irreversible, we warmed 
the crystal above Td and re-cooled to T = 60 K in zero 
field. We then swept the field up and down several times, 
reaching successively larger values each time. We found 
that any increase in susceptibility, no matter how small, 
was irreversible (see Fig. I12f b % l). For H || [110] we found 
the susceptibility to be completely independent of field 
at T = 60 K. The data in Fig. 1121^ have been corrected 
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FIG. 11: Field-dependence of the susceptibility for field ap- 
plied at constant T = 2K. (a) H || [001], (b) H || [110]. The 
small upturns in the susceptibility at low fields in both (a) and 
(b) are probably due to underestimation of the applied field 
strength, due to magnetisation of the sample surroundings. 



FIG. 13: Temperature-dependence of M/H at H = 0.3 T for 
H || [001] before and after application and removal of a 7T 
field, (a) M/H vs T. The two traces coincide above Trj. (b) 
{M/H)' 1 vs T. The anomaly at T D is more pronounced after 
application and removal of the 7T field. 



for a remanent field of = 0.008 T, which was prob- 
ably caused by magnetisation of the sample surround- 
ings or trapped flux in the superconducting magnet. Hr 
was determined by finding the y-intercept of a straight 
line fitted to a plot of magnetisation M versus field H 
(M = jj-{H + -Hr), where .Hr is the remanent field). 
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FIG. 12: Field-dependence of the susceptibility for H || [001], 
field applied at constant T = 60 K. The data in (a) has been 
corrected for a small remanent field of Hb. = 0.008 T, prob- 
ably caused by magnetisation of the sample surroundings or 
trapped flux in the magnet, (b) shows that any increase in 
susceptibility, no matter how small, is irreversible. 



The irreversible increase in susceptibility for H || [001] 
persists up to Td. This is shown in Fig. ll3f a). where the 
temperature dependence of the susceptibility, measured 
at H = 0.3 T is plotted before and after application and 
removal of a 7T field. The two traces coincide for T > 
Tjy. Figure ITST b - ) shows the inverse susceptibility (taken 
from the same data set as Fig. I13f a)). The anomaly at 
Td is more pronounced after application and removal of 
the 7T field. 



C. Summary 

The main results of the susceptibility studies are as fol- 
lows. The anomaly at Td is more pronounced when the 
field is applied along [110] than when it is applied along 
[001]. For H || [110] the susceptibility increases reversibly 
with field below Tn, but is independent of field above 



Tn. For H || [001] the susceptibility increases rapidly be- 
low Tn from i? = 0Ttoi/ = lT and more gradually 
from H = IT to H = 7T. The increase from 0-1 T is 
quasi-reversible and occurs only below Tn, whereas the 
increase from 1-7 T is irreversible, and also occurs above 
Tn. The field strength at which the latter occurs depends 
on whether the field is applied below or above Tn. For 
fields applied at constant T < Tn a field of 7 T is required 
to achieve the full increase in susceptibility, whereas for 
fields applied above Tn, only 5.5 T is required. Any in- 
crease in susceptibility above Tn , no matter how small, is 
irreversible. Above Td the susceptibility is independent 
of field. 



V. INTERPRETATION OF SUSCEPTIBILITY 
DATA 

We now attempt to interpret the susceptibility data 
using our hypothesis that the formation of magnetic do- 
mains is strongly influenced by the distribution of struc- 
tural domains. 

The most likely cause of the reversible increase in sus- 
ceptibility with field for H || [110] is the proposed re- 
versible spin reorientation in the type-I component of the 
magnetic structure. 

For H | [001] we propose that the rapid quasi- 
reversible increase in susceptibility between H — T and 
H = 1 T is due to initial favouring of magnetic domains 
whose spins are perpendicular to the field (D b a , Dg, 
and D b c ). This would give rise to a small decrease in 
the magnetic structure factor of the (100) Bragg peak 
and a small increase in the magnetic structure factor of 
the (110) Bragg peak, in reasonable agreement with the 
small changes in Bragg intensities observed at H = 1 T 
(see Fig. Between H = IT and H = 7T we pro- 
pose that the magnetic field becomes large enough to 
influence the alignment of the structural domains. The 
resulting structural domains cause those type-I magnetic 
domains whose ordering vectors are perpendicular to H 
to be favoured. This leaves only domains D h a and 
causing the magnetic structure factors, and hence the in- 
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tensities of the Bragg reflections, to change. When the 
field is decreased below 1 T, we believe that domains D c a 
and Dl return. Combining the structure factors of these 
domains with those of D h a and makes no difference to 
the average structure factors, in agreement with the ob- 
served absence of change in Bragg intensities on removal 
of the field. However, the return of these domains with 
spins parallel to the applied field would be expected to 
cause a decrease in susceptibility, as observed. 

The above discussion shows that the hypotheses out- 
lined in Scction lTTTl Drovide some success in explaining the 
susceptibility data. However, some observations remain 
unexplained. For instance, we observe an irreversible in- 
crease in susceptibility between H = 1 T and H = 7 T 
for H || [001] both below and above Tn. We assume that 
the increase above Tn is due to the effect of the applied 
field on the ordering of the Pr orbitals. However, no ir- 
reversible increase is observed either above or below Tn 
for H || [110], and we do not understand why this should 
be. We also do not understand why a change in domain 
populations between H = T and H = IT should cause 
a change in susceptibility for H || [001] when the irre- 
versible change in domain populations after application 
of a large field along [110] causes no irreversible change 
in the susceptibility below Tn- 

Despite the above problems, the interpretation of the 
susceptibility data is in broad agreement with that of 
the neutron diffraction data. We also note that the sus- 
ceptibility measurements at low fields indicate that the 
[001] direction is the easy direction of magnetisation, in 
support of a single-q type-I component of the magnetic 
structure. 



VI. CONCLUSION 

We have presented (i) neutron diffraction studies of the 
antiferromagnetic structure of single crystal Pr02 as a 
function of applied field (ii) susceptibility measurements 



on single crystal PrC>2 as a function of applied field and 
temperature. Both studies were carried out for fields 
applied along two different symmetry directions: [001] 
and [110]. 

The neutron diffraction studies revealed irreversible 
changes in the magnetic Bragg intensities after applica- 
tion and removal of a large field. We attributed these 
changes to magnetic domain alignment, caused by a field- 
induced change in the underlying distribution of struc- 
tural domains. For H || [011] we also observed certain 
reflections to disappear at high field and reappear when 
the field was removed. We attributed this observation 
to a spin reorientation of the type-I component of the 
magnetic structure. 

The susceptibility studies revealed a reversible increase 
with applied field below Tn along both symmetry direc- 
tions, as well as an irreversible increase at higher fields 
both above and below Tn for H || [001]. The irreversible 
increase disappeared at Td, the temperature at which 
a cooperative Jahn- Teller distortion of the oxygen sub- 
lattice occurs. The changes in susceptibility with applied 
field are in broad agreement with the proposed alignment 
of structural and magnetic domains, although some fea- 
tures remain poorly understood. 

We conclude that our neutron diffraction and suscep- 
tibility studies have revealed striking hysteresis effects 
on the application and removal of a large magnetic field 
below Td- We propose that the field interacts strongly 
with the lattice to influence the alignment of structural 
domains below Td, and that this in turn influences the 
alignment of the magnetic domains below Tn- 
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